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ABSTRACT 


An  investigation  of  using  the  response  of  a  circuit  at  selected  test 
frequencies  to  isolate  faulty  circuit  components  is  made.  A  procedure 
using  a  sensitivity  approach  for  intelligent  selection  of  test  frequencies 
is  developed.  The  developed  procedure  is  tested  and  the  results  are 
compared  with  results  using  conventional  procedures.  The  effect  of 
random,  within  tolerance  variations  of  nor.faulty  components  on  the 
results  is  studied  for  >oth  conventional  and  developed  procedures.. 
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t .  TNTRnrmrTTnN 


Maintenance  of  electronic  equipment  is  a  problem  that  increases  as 
equipment  becomes  more  complex.  With  ever  increasing  costs  (of  equip¬ 
ment)  it  also  becomes  imperative  to  decrease  equipment  downtime  as 
much  as  possible.  To  this  end  a  fast,  reliable  method  of  locating 
circuit  failures  must  be  used. 

It  has  been  customary  to  isolate  failures  by  testing  individual  compo¬ 
nents.  This  is  usually  done,  by  isolating  the  fault  to  a  particular  stage, 
followed  by  subassembly  isolation  and  finally  component  isolation  by 
individual  testing. 

With  the  availability  and  flexibility  of  digital  computers,  a  pos¬ 
sible  solution  to  the  problem  of  excessive  downtime  is  to  use  a  computer 
to  help  isolate  the  fault.  The  first  step  in  accomplishing  this  is  the 
development  of  a  testing  and  isolation  procedure  which  could  be  program¬ 
med  for  the  digital  computer. 

A  recently  presented  procedure  has  been  examined  by  a  number  of 
individuals  [1],  [2],  [3],  and  [4].  The  procedure  was  described  by  Seshu 
and  Waxman  [1]  with  other  investigators  testing  and  making  applications  of 
the  procedure  with  few  modifications. 

The  basic  procedure  is  to  examine  the  response  of  a  circuit  under  test 
at  selected  frequencies.  Comparisons  of  the  results  with  a  list  of  “library 
results"  can  pinpoint  a  faulty  component. 

Most  investigators  have  used  test  frequencies  near  the  poles  and  zeros 
of  the  function  under  consideration.  Other  frequencies  have  been  chosen 
between  the  pole  and  zero  frequencies  [2]  and  [4],  In  this  study  an 
alternate  method  of  test  frequency  selection  was  developed  and  tested  and 
results  were  compared  with  results  using  other  methods. 
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The  procedures  for  signature  comparison  used  by  other  investigators 
have  not  been  explained  in  detail  in  most  reports.  However,  it  appears 
that  a  nearest  neighbor  approach  is  the  most  conation  technique  used  to 
date.  Two  alternate  methods  were  examined  and  are  discussed  with 
conclusions  as  to  their  practicality. 

Previous  methods  of  fault  determination  after  signature  comparison 
have  also  not  been  documented.  Two  possible  methods  were  tried.  The 
first  method  used  computer  selection  for  each  network  function  and  manual 
correlation  for  fault  isolation.  The  second  method  was  completely  program¬ 
med  for  digital  computation. 

As  far  as  can  be  determined  from  available  reports  no  investigator 
approached  the  problem  of  error  due  to  in-tolerance  variations  of  com¬ 
ponents.  This  study  included  an  investigation  of  that  problem  with  both 
the  conventional  and  the  proposed  procedures. 
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I I .  PROCEDURAL  BACKGROUND 

A.  FUNDAMENTAL  CONCEPT 

All  components  of  a  given  network  play  a  part  in  the  formation  of  any 
of  the  network  functions.  The  sensitivity  of  any  given  network  function 
to  a  particular  component  is  generally  unique  to  that  component.  That  is, 
in  a  network  with  K  components  a  network  function  will  have  K  sensitivity 
functions,  one  for  each  component.  A  change  in  the  value  of  a  particular 
component  may  cause  a  decrease  in  the  magnitude  of  a  network  function  at 
one  frequency  and  an  increase  in  the  magnitude  at  a  different  frequency. 
There  may  in  fact  be  frequencies  at  which  the  function  sensitivity  is 
near  or  equal  to  zero  for  a  particular  component,  excluding  catastrophic 
failures. 

Specific  knowledge  about  this  behavior  can  be  used  to  detect  and 
Isolate  component  value  changes.  If,  for  example,  it  is  known  that  the 
Input  impedance  is  insensitive  at  low  frequencies  to  changes  in  a  specific 
component,  one  would  assume  that  component  was  not  faulty  if  the  input 
impedance  exhibits  a  significant  change  at  Igw  test  frequencies. 

B.  POLE-SHIFT  TECHNIQUE 

The  component  variation  effect  on  frequency  response  can  be  observed 
with  the  pole-shift  technique.  A  network  function  can  be  represented  by 
a  constant  and  a  unique  set  of  poles  and  zeros;  i.e., 

F(S)  -  k  fr-yfr-y  •••  <s-y 

K  {s-p-,j(s-p2)  ...  (s-pn) 

where 

F($)  *  any  network  function 


K  *  constant 


-th 

Z|  *  1’"  zero  frequency 

.  fh 

Pj  *  y"  pole  frequency. 

Every  component  plays  a  part  in  pole  and  zero  locations.  However*  a 
few  components  may  dominate  the  location  of  a  given  pole  or  zero  frequency. 
Thus,  if  one  could  detect  the  shift  of  the  poles  and  zeros  of  a  given 
function,  it  may  be  possible  to  determine  component  variations.  A  simple 
example  can  best  illustrate  this  point. 

Example  1  -  Consider  the  network  shown  in  Figure  1. 
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Figure  1.  Example  1  -  Low  Pass  Filter 


The  network  function 


(s)  is  given  by: 


V1  _  s2  +  11s  +  11 
_  *  - - 

*1  s*  +  103  +  1 


s  +  9.887) (s  +  1.113) 
s  +  9.899) (s  +  0.1 01}  * 


The  pole  and  zero  shifts  were  determined  by  varying  each  component 
value  by  ten  percent  and  observing  the  corresponding  critical  frequency 
shifts.  The  results: 
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Location  (radians/sec) 


K  - a.  AL  J  Jf>JL 

rercent  orurt 


Condition 

h 

P2 

21 

z2 

Pi 

P2 

21 

Nominal 

.101 

9.899 

1.113 

9.887 

0 

0 

0 

0 

R.j  variation 

.101 

9.899 

1.020 

9.889 

0 

0 

.62 

.02 

Cj  variation 

.092 

9.908 

1.010 

9.899 

8.9 

.09 

9.25 

.12 

L^  variation 

.101 

8.990 

1.114 

8.977 

0 

9.2 

.08 

9.2 

R2  variation 

.092 

10.908 

1.101 

10.899 

8.9 

10.2 

1.07 

10.2 

TABLE  I 1-1.  Percent  Pole  and  Zero  Shifts 


It  Is  immediately  obvious  that  a  change  in  R^,  or  C-j  will  cause  a  greater 
shift  in  z-j  than  will  a  corresponding  change  in  Lj  or  R^.  Also,  p2  and  z2 
are  most  affected  by  L-j  and  R2  and  p-j  is  affected  only  by  and  R^ 

C.  CHANGES  IN  RESPONSE  FROM  POLE  SHIFT 

In  practice  the  exact  location  of  the  critical  frequencies  can  usually 
not  be  determined  directly  from  the  circuit.  However,  this  does  not  rule 
out  the  pole  shift  technique.  Any  shift  in  critical  frequency  causes  a 
change  in  the  magnitude  of  the  given  function  at  all  frequencies  greater 
than  the  shifted  one.  This  can  be  seen  from  the  simple  Bode  representation 
of  a  pole: 
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If  f  Is  changed  to  f'  (indicated  by  dashed  line),  the  function 
magnitude,  at  any  frequency,  fQ,  greater  than  fc,  is  changed  (increased 
in  this  case).  The  magnitude  is  not  affected  at  frequencies  less  than 

v 

Now,  from  the  network  in  Figure  1  we  obtain: 

V2  c  10  =  10 _ 

^  s2  +  11s  +  11  (s  +  9.S87)(s  +  1.113) 


Figure  3a  is  the  Bode  plot  for  the  function  ^2^1'  ^9ure  and 
3c  show  the  effect  of  increasing  R-j  and  of  increasing  ,  respectively. 
The  nominal  response  is  shown  as  a  dashed  line. 
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It  can  be  seen  that  an  increase  in  R.,  causes  a  decrease  in  the 
magnitude  of  the  response  between  and  Pg.  A  lesser  decrease  Is  seen 
at  frequencies  greater  than  pg.  The  response  between  p^  and  pg  is 
increased  for  an  increase  in  Lj  while  the  response  above  Pg  is  decreased. 
Thus,  if  one  examines  the  response  at  a  frequency  between  p^  and  p^  and 
at  a  frequency  greater  than  Pg  it  will  be  a  reasonably  simple  task  to 
distinguish  between  an  increase  in  R-j  or  l-j.  This  example  is  an  over¬ 
simplification  of  the  procedure,  but  the  application  of  which  can  be 
appreciated. 
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III.  GENERAL  PROCEDURE 

A.  SASIC  PROCEDURE 

The  basic  procedures  outlined  by  Seshu  and  Waxman  [1]  was  followed 
with  modifications  as  explained  herein.  Two  network  functions  were  used 
concurrently.  Work  done  by  Maenpaa,  Stehman,  and  Stahl  [5]  indicated 
that  the  redundancy  from  using  two  functions  was  more  than  adequately 
compensated  for  by  the  increased  fault  detection  capability. 

The  symbolic  network  functions  were  calculated  by  hand.  It  should 
be  pointed  out  that  computer  programs  exist  which  calculate  symbolic  net¬ 
work  functions.  These  programs  should  be  used  for  more  complex  circuits. 

The  functions  were  numerically  evaluated  and  factored  to  determine 
the  poles  and  zeros.  Test  frequencies  were  then  chosen.  Previous 
procedures  generally  selected  test  frequencies  based  on  the  pole  and 
zero  locations.  As  a  test  of  the  procedure  several  test  frequencies  were 
chosen  between  the  critical  frequencies  as  well  as  0  hertz  and  one 
frequency  above  the  highest  critical  frequency.  The  actual  selection 
was  programmed  for  computer  selection  and  is  described  in  Section  VII. 

It  was  felt  that  the  procedure  of  selecting  test  frequencies  based 
on  poles  and  zeros  was,  although  intuitively  satisfying,  a  somewhat 
illogical  procedure.  The  concern  was  not  with  the  function  itself  but 
with  the  behavior  of  that  function  with  a  parametric  change.  Thus, 
sensitivity  functions  were  calculated  for  each  element.  Based  on  these 
functions,  six  test  frequencies  were  chosen.  Th.s  procedure  is  discussed 
In  Section  IV. 

After  selection  of  test  frequencies,  a  worst  case  response,  for 
component  values  within  10  percent  tolerance  limits,  was  computed  for 
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each  test  frequency.  This  defined  a  nominal  range  of  performance.  This 
is  not  necessarily  the  best  method  for  definition  of  nominal  range.  This 
range  might  be  specified  for  a  particular  circuit.  However,  in  the 
absence  of  such  specifications  the  above  procedure  was  followed. 

Each  component  was  then  allowed  to  take  on  several  discrete  values 
out  of  tolerance  limits.  All  other  components  were  held  fixed  at  their 
nominal  value.  The  response  at  each  test  frequency  was  then  calculated 
and  catalogued  for  the  entire  collection  of  sets  of  component  values. 

The  catalogued  response  values  were  then  quantized  into  nine  levels 
and  labeled  one  through  nine.  Quantization  level  five  was  reserved  for 
response  within  nominal  limits.  The  setting  of  the  levels  is  explained 
in  Part  B.  The  quantized  responses  for  a  given  set  of  component  values 
is  defined  as  the  signature  for  the  given  circuit  configuration. 

The  signatures  were  computed  for  the  following  cases: 

value  =  zero 
♦  2W 
+  50X 

100  times  nominal. 

The  entire  set  of  signatures  for  a  given  function  is  defined  as  the 
library  for  that  function.  There  is  no  theoretical  limit  on  the  size 
of  the  library.  However,  the  inclusion  of  many  more  circuit  configurations 
would  probably  only  complicate  the  matching  procedure  and  not  Increase 
detection  capability. 

Using  the  primary  test  circuit  the  input  Impedance  and  voltage  gain 
at  each  test  frequency  was  measured.  Then,  after  quantization,  the  signa¬ 
tures  measured  were  compared  with  the  library  signatures  and  a  "closest 
match"  was  made.  Thus,  the  faulty  component  was  identified. 

In  order  to  test  the  procedure,  several  error  conditions  were  set  and 
the  response  calculated.  In  order  to  more  closely  approximate  actual 
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B.  SETTING  OF  QUANTIZATION  LEVELS 

With  the  exception  of  selection  of  test  frequencies  the  most  critical 
step  in  the  procedure  was  the  setting  of  quantization  levels.  Two  methods 
were  used  successfully.  One  was  programed  as  part  of  the  main  analysis 
program.  The  second  method,  although  not  programed  is  possibly  adaptable 
to  computer  selection.  Both  methods  have  nine  levels  of  response  with 
level  five  reserved  for  the  nominal  range. 

The  nominal  range  was  selected  by  taking  the  range  between  the 
maximum  positive  and  maximum  negative  deviation  from  normal  for  any 
element  value  at  +  10%  from  normal.  The  programmed  method  then  set  the 
other  eight  levels  as  shown  in  Table  III-l  in  which  F  is  a  multiplying 
factor  used  to  force  the  possible  response  range  to  cover  all  quantum 
levels  with  levels  ont  and  nine  assigned  only  for  extreme  cases. 


Quantum  Level 


Normalized  Response 


1 

2 

3 

4 

5 

6 

7 

8 
9 


(-,  1  +  13  FN] 

(1  +  13  FN,  1  +  7  FN] 
(1  +  7  FN,  1  +  4  FN] 
(1+4  FN.  1  +  N] 

(1  +  H,  1  -  M) 

[1  -  N,  1  -  4'FM) 

[1  -  4  FM,  1  -  7  FM) 
[1-7  FH,  1-13  FM) 
[1-13  FM,  0) 


TABLE  III-l.  Quantization  levels. 
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deviation  levels,  1  and  9,  were  selected  based  on  the  maximum  deviation 
from  normal  under  any  condition.  The  mid  levels  were  then  selected  so 
that  the  catalogued  or  library  signatures  would  cover  all  quantum  levels 
as  uniformly  as  possible. 

The  manual  approach  is  acceptable  for  small  circuits  but  gets 
unmanageable  quickly  as  circuit  size  increases.  Parts  of  the  procedure 
could  be  easily  programmed.  However,  a  suitable  algorithm  would  have  to 
be  developed  prior  to  complete  programming. 


C.  TWO  COMPONENT  VARIATIONS 

It  cannot  be  assumed  that  only  one  component  at  a  time  will  undergo 
changes  in  value.  In  order  to  investigate  the  possibility  of  two  compo¬ 
nent  values  changing  simultaneously,  two  possibilities  were  considered. 

If  one  comoonent  masks  function  changes  due  to  variations  in  another 
component  (i.e.,  no  change  in  function  magnitude),  the  entire  procedure 
will  not  work.  However,  if  one  is  primarily  concerned  with  circuit 
performance,  the  failure  to  detect  component  variations  is  of  no  concern 
If  there  is  no  degradation  in  circuit  performance. 

The  second  situation  considered  was  one  in  which  there  was  an  obvious 
change  in  circuit  response  due  to  two  components  varying.  The  original 
procedure  made  allowances  for  varying  only  one  component  at  a  time.  How¬ 
ever,  a  slight  modification  allowed  for  the  generation  of  a  set  of  library 
signatures  fGr  two  component  variations.  The  set  is  presented  in  the 
computer  output  section. 

It  soon  became  apparent  that  in  order  to  use  the  two  component 
variation  library,  one  needed  apriori  knowledge  of  the  number  of  faulty 
components.  Since  this  is  not  usually  available,  the  usefulness  of  two 
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component  signatures  is  doubtful.  Additionally,  an  examination  of  the 
I  lurarjr*  iiynawi cs  TcVcalcd  that  One  of  the  tWO  COiTipOnentS  Varied  USUally 
dominates  the  response.  Thus,  if  two  components  have  failed  the  single 
component  library  would  probably  isolate  the  dominant  component.  After 
replacement  of  that  component  further  testing  could  reveal  the  second 
faulty  component.  Consequently  no  further  investigation  of  multi¬ 


component  variation  was  undertaken. 
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IV.  PRIMARY  TEST  CIRCUIT 


The  primary  circuit  used  throughout  the  investigation  was  a  six- 
element  low-pass  filter.  (See  Figure  4.)  The  unsealed  frequency 
response,  represented  by  the  voltage  ratio  transfer  function,  V2/V-j* 
is  shown  in  Figure  5  and  the  unsealed  input  impedance  as  a  function 
of  frequency  is  shown  in  Figure  6.  On  Figures  5  and  6  the  poles  are 
represented  by  X  and  the  zeros  by  0. 


The  component  values  were  magnitude  scaled  by  100  and  frequency 

c 

scaled  by  10  .  This  yielded  component  values  of  approximately  the 
same  magnitude.  Ideally,  the  scaled  cutoff  frequency  should  be  1  Hz. 
However,  with  this  scaling  f  *  .02  Hz  and  no  difficulty  was  encounter- 

V» 

ed. 

The  scaled  network  functions  of  interest  are: 

V2  ,  -1634 _  ‘ 

^7  s4+l 7 . 67s3+6 . 462s2+2 . 601 s+ . 21 24 

.  _ _ .1634  _ 

(s+1 7 . 31 ) (s+ . 09$7 ) ( s+i . 325+j 3.26) (s+1 .325-J3.2B) 
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wit  gegpfwg  ggg aggggwgp  '}  m  i  ■  jh  lj  mm  ,J'-ujj  g»  m  m*  w  *p*  Ju^F^-g  w 


V1  ,  s4-H  7 . 67s3t-6 . 462s2+2 . 601  s+ .  21 24 
*7  5.882s3+.147s2+.817s+.0] 634 


(s+17.31 )(s+.0997) 

(s+1 .325+j3.25){s+l .325-j3.25) 

(s+.O2)(s+.0C25+j. 

573)(s+.0d25-j.3?3)“ 

The  poles  and  zeros  in  Hz  are  listed  in  Table  IV-1  for  convenience. 


V*i 

V*i 

Zeros* 

(Hz) 

Poles  (Hz) 

.0518 

(Second  Order) 

.0594  (Second  Order) 

.0159 

.00319 

2.76 

*  Zeros  of  V^/I^  are  the  poles  of  V2/V^. 


TABLE  IV-1.  Pole  and  Zero  Magnitudes. 


V.  SELECTION  OF  TEST  FREQUENCIES 


A.  CONVENTIONAL  METHOD 


The  conventional  procedure  was  followed  with  regard  to  selection  of 
test  frequencies  in  order  to  test  the  method.  It  soon  became  apparent, 
that  the  number  of  test  frequencies  was  too  large  to  handle  with  ease. 
Thus,  only  frequencies  between  the  pole  and  zero  frequencies  were 
selected  as  well  as  0  Hertz  and  one  above  the  maximum  critical  frequency. 

The  computer  program  written  for  implementation  of  the  procedure  (see 
Program  1)  did  the  frequency  selection  automatically.  The  test  fre¬ 


quencies  were  selected  by  taking  4,, 

.  fl  1  -  f 


f  =  +  c  , 

Ttest  Tc  -  T 


ci  .  .ci  _  -*th 


for  every  f'  where  f  =  i  critical  frequency.  If 


\  f !+1  -  fi  ' 
fi  +  c  c 

xi+2  ,i+l 
*i+l  Tc  "  fc 

Lc  2  j 

* 

fc  ~  2 

one  of  the  frequencies  was  eliminated.  Test  frequencies  were  computed 
for  the  input  impedaw-e  because  the  critical  frequencies  for  V2/V,  are 
Included  in  the  critical  frequencies  for  V-j/Ij.  One  test  frequency  was 
eliminated  due  to  the  above  inequality.  This  resulted  in  the  selection 
of  the  following  seven  test  frequencies  for  a  fourth  order  filter: 


0.0 

Radians 

.001 

Radians 

.00498 

Radians 

.196 

Radians 

.5488 

Radians 

8.49 

Radians 

25.8 

Radians 

TABLE  V-l .  Conventional  Method  Test  Frequencies 
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B.  SENSITIVITY  APPROACH  TO  FREQUENCY  SELECTION 

Inspection  of  the  signatures  generated  using  the  conventional  method 
showed  that  redundancy  existed  to  the  extent  that  the  method  was  some¬ 
what  inefficient.  Two  test  frequencies,  D.C.  and  the  lowest  nonzero  test 
frequency,  gave  nearly  identical  results.  This  indicated  one  of  the  two 
frequencies  was  probably  not.  required.  It  was  also  felt  that  unless 
further  investigation  was  undertaken  it  would  be  possible  to  ignore  a 
range  of  frequencies  which  could  be  extremely  helpful. 

* 

This  investigation  took  the  form  of  calculating  and  plotting  functions 
which  defined  the  sensitivity  of  Vg/V-j  and  I-j/V-j  to  variations  in  each 
component.  These  plots  are  presented  as  Figures  7  through  18.  Figure  19 
is  a  superposition  of  the  plots  of  the  sensitivity  function  for  I^/V^  and 
Figure  20  is  the  same  for  Vg/V-j .  Ideally  V^/I^  should  have  been  used. 
However,  the  computer  program  used  in  this  portion  of  the  study  did  not 
allow  for  using  Vj/I-j.  It  vas  felt  that  no  loss  of  information  would 
result. 

The  superposition  plots  Figures  19  and  20  are  very  revealing.  The 
first  observation  is  that  V^/V-j  is  not  sensitive  to  changes  in  R-j  and 
at  high  frequencies  nor  to  C's  or  L's  at  low  frequencies.  While  this 
information  is  not  new  the  sensitivity  functions  assist  in  determining 
when  C's  and  L's  have  little  effect.  This  can  be  used  to  eliminate  the 
redundancy  observed  using  the  conventional  method  of  frequency  selection. 

The  mid -frequency  range  is  the  most  interesting  range  for  the  selection 
of  test  frequencies.  For  example,  near  .02  Hz  Vg/V-j  is  very  sensitive  to 
changes  i.i  R^  and  less  sensitive  to  changes  in  Lg,  even  less  sensitive 


*  The  program  NASAP-70  was  used  to  calculate  the  various  sensitivity 
functions. 
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to  changes  In  R2  and  Cj  and  virtually  insensitive  to  changes  in  Lj. 

There  is  a  peaking  of  the  sensitivity  functions  for  Cj,  C2,  and  L,  near 
.06  Hz.  At  this  same  frequency  the  sensitivity  functions  for  Rj,  R2» 
and  L|  are  near  zero. 

The  I/.j  and  V-j  plots  show  that  at  near  .3  Hz  only  the  sensitivity 
function  for  Rj  is  high.  The  functions  for  Lj  and  Cj  are  low  and  the 
other  functions  are  zero.  At  3  Hz  Ij/Vj  is  about  equally  sensitive  to 
changes  in  Lj  and  Rj  while  the  other  functions  are  very  close  to  zero. 

It  can  also  be  seen  that  the  frequencies  discussed  in  connection  with 
V2/Vj  are  reasonably  good  choices  for  Ij/Vj. 

Based  on  the  above  observations  the  following  -lx  frequencies  were 
selected  as  test  frequencies: 

Frequency 

.00048  Hz 
.027  Hz 
.06  Hz 
.299  Hz 
2.99  Hz 
49.97  Hz 

TABLE  V-2.  Sensitivity  Function  Test  Frequencies. 

After  some  preliminary  work  it  was  observed  that  the  response  at  .06 
Hz  was  not  well  behaved.  Since  four  sensitivity  functions  peaked  for 
Ij/Vj 9  random  variation  within  tolerance  limits  of  these  components  were 
being  observed  as  significant  changes  in  response.  In  order  to  overcome 
this  problem,  test  frequency  three  was  changed  to  .055  Hz.  The  new  fre¬ 
quency  was  r.ff  the  peak  enough  to  be  useable  and  still  yield  the  desired 
information. 


Number 

1 

2 

3 

4 

5 

6 


4-|iF.(w)| 


3. 


2. 


1. 


0.1 _ 

.001 


.01  .1 
Figure  9.  Sensitivity  of 
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.001  .01  .1  1.  10.  w 

Figure  15.  Sensitivity  of  Ij/V^  to  Cl. 
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Figure  16.  Sensitivity  of  I^/V^  to  R, 


•  Figure  19.  Ij/V-j  Composite  Sensitivity  Functions. 
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C.  COMPARISON  OF  METHODS 


The  comparison  of  the  results  using  the  two  methods  of  frequency 
selection  is  included  in  Section  VI.  Some  comments  as  to  applicability 
should  be  made  at  this  point. 

If  one  computes  sensitivity  functions  by  using  the  derivative  of 
the  network  function  with  respect  to  the  element  under  consideration  it 
becomes 

fs(s)  =  -ai^ 


«w 

_ k  _ k_ 

~  Ms)]2 


where: 


fs(s)  *  sensitivity  function 

d/de^  =  derivative  with  respect  to  element  k 


»  ratio  of  polynomials  representing  the  network  function. 


Thus,  the  poles  of  the  sensitivity  function  are  the  same  as  the  poles  of 
the  original  function.  Consequently,  if  one  selects  test  frequencies  on 
the  basis  of  function  poles  the  selection  is  linked  to  the  sensitivity 
functions.  Generally,  however,  the  zeros  of  the  sensitivity  function  are 
unique  and  are  thus  ignored  by  any  selection  procedure  based  on  the  poles 
and  zeros  of  the  original  function. 

From  the  preceding  paragraph,  some  of  the  advantages  of  using  sensi¬ 
tivity  functions  for  test  frequency  selection  become  apparent.  However, 
that  method  has  the  disadvantage  of  requiring  the  calculation  of  sensi¬ 
tivity  functions  for  each  element. 
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With  the  small  circuit  used  as  an  example  the  value  of  sensitivity 
functions  are  demonstrated.  With  more  complex  circuits  it  appears  more 
probable  that  important  test  frequencies  might  be  overlooked  if  the  sen 
sitivity  functions  are  not  used. 


VI.  RESULTS  WITH  PRIMARY  CIRCUIT 

A.  RESULTS  USING  BASIC  PROCEDURE 

Using  the  procedures  outlined  in  Section  III-A  and  V-A  a  set  of 
signatures  was  generated  for  the  six  selected  circuit  conditions. 

These  signatures  are  presented  in  Figure  25.  A  second  set  of  signatures 
was  then  generated  for  the  same  circuit  conditions  with  random  variation 
allowed  for  the  circuit  elements.  The  second  set  of  signatures  was  input 
to  the  selection  program  to  determine  what  errors  resulted  from  the  random 
variations. 

Table  VI -1  lists  the  input  condition  and  the  program  output.  The 
numerical  result  was  determined  as  follows: 

1  =  Selection  was  correct  condition 

2  =  Selection  was  correct  element 

3  =  One  of  multiple  selections  was  the  correct  condition 

4  *  One  of  multiple  selections  was  correct  element 

5  =  Incorrect  selection. 

By  correlating  the  results  from  Vg/V-  and  V-j/Ij  a  “hit",  "miss”,  or 
"possible  hit"  rating  was  assigned  for  each  test  signature  as  follows. 

A  hit  was  assigned  if  both  V^V-j  and  Vj/Ij  yielded  a  1  or  2  or  if  one 
was  1  or  2  and  the  other  was  3  or  4.  A  3-5,  4-5,  or  5-5  combination  was 
called  a  miss.  All  other  combinations  were.called  possible  hits.  Fcr 
some  cases  deviations  were  warranted.  For  example,  a  3-4  selection  could 
have  been  called  a  miss  if  a  component  other  than  the  given  component 
was  clearly  indicated  as  being  faulty. 

The  result  was  21  hits,  9  possible  hits  and  6  misses  for  58%  hits  and 
83%  hits  or  possible  hits.  Since  only  the  circuit  conditions  used  for 
signature  generation  were  used  to  generate  the  test  signatures  100% 
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Circuit 

Program 

Numerical 

Program 

Numerical 

Condition 

Condition 

Result 

Condition 

Result 

R1:  Short 

Rl:  Short 

1 

Rl:  Short 

1 

-50* 

Rl:  -50* 

1 

L2:  Open 

5 

-20* 

Rl:  Open 

2 

L2:  +50* 

5 

+20* 

L2:  +20* 

5 

Rl;  +20* 

1 

+50* 

Rl:  +20* 

2 

Rl:  +20* 

2 

Open 

Rl:  Open 

1 

L2:  Short 

Rl:  Open 

3 

LI:  Short 

LI:  Short 

1 

LI:  -50* 

2 

-50* 

LI:  -50* 

1 

LI:  -50* 

1 

-20* 

LI:  -20* 

R2:  +20* 

3 

LI:  -20* 

1 

+2C* 

LI:  +20* 

1 

LI:  -20* 

2 

+50* 

LI :  +20* 

2 

LI:  -20* 

2  ' 

Open 

LI:  Open 

1 

R2:  +20* 

5 

Cl:  Open 

Cl:  Open 

1 

LI:  +50* 

5 

-50* 

Cl:  +20* 

C2:  -50* 

4 

Cl:  -50* 

1 

-20* 

C2:  +20* 

5 

Rl:  -20* 

5 

+20* 

LI :  +20* 

C2:  +20* 

.  5 

LI:  -50* 

5 

+50* 

C2:  +20* 

5 

Cl:  +20* 

C2:  +20* 

3 

Short 

Cl:  Short 

1 

Cl:  Short 

1 

R2:  Short 

R2:  -50* 

2 

R2:  Short 

t 

-50* 

R2:  -50* 

1 

R2:  -50* 

R2:  -20* 

2 

-20* 

R2:  -20* 

L2:  +20* 

C2:  +20* 

3 

Rl:  -20* 

5 

+20* 

C2:  +20* 

5 

Cl:  +20* 

C2:  +20* 

5 

+50* 

R2:  +20* 

2 

R2:  +20* 

2 

Open 

R2:  Open 

1 

Cl :  Open 

5 

TABLE  VI -1.  Basic  Procedure  Results. 
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L2:  Short 

L2:  Short 

1 

L2:  Short 
Rl:  Open 

-50% 

L2:  -50% 

1 

L2:  -50% 

-20% 

L2:  -50% 

LI:  +20% 

2 

L2:  -20% 
L2:  +20% 

+20% 

L2:  +20% 

C2:  +20% 

3 

L2:  Open 

+50% 

L2:  +20% 

2 

Open 

L2:  -20% 

2 

L2:  Open 

C2:  Open 

Cl:  +20% 

C2:  -50% 

4 

R2:  +20% 

-50% 

Cl:  -20% 

5 

Cl:  +50% 
Cl:  -20% 
C2:  -20% 

-20% 

Nominal 

5 

L2:  +20% 

+20% 

C2:  +20% 

1 

C2:  +20% 

+50% 

Nominal 

5 

C2:  +20% 
Cl:  +20% 

Shore 

C2:  Short 

R2:  Short 

3 

C2:  Short 

TABLE  VI -1  (Continued) 


hits  or  possible  hits  could  have  been  accomplished.  Thus,  using  the 
conventional  procedure,  within  tolerance  variations  introduced  172  to 
422  errors. 

B.  RESULTS  USING  SENSITIVITY  FUNCTIONS 

Using  the  frequencies  determined  from  the  sensitivity  functions  a 
set  of  signatures  was  generated  for  the  six  circuit  conditions.  In  order 
to  check  the  effect  of  variations  within  tolerance  for  this  method  the 
same  procedure  as  used  in  VI-A  was  followed.  The  result  was  50  hits,  16 
possible  hits  and  6  misses.  The  6  misses  included  3  due  to  the  test 
signature  showing  normal  response  from  the  circuit.  These  results  showed 
an  error  rate  of  112  -  312. 

Eight  additional  circuit  conditions  were  simulated  for  each  element 
to  see  what  could  be  expected  from  a  less  artificial  situation.  Table 
VI-2  lists  all  the  conditions  set  and  the  numerical  result.  The  numerical 
result  was  assigned  as  follows: 

1  -  closest  library  signature  selected 

2  -  right  element  selected 

3  -  one  of  multiple  signature  closest  to  given  condition 

4  -  one  of  multiple  signatures  right  element 

5  -  incorrect  selection. 

Hits,  misses  and  possible  hits  were  also  assigned  as  before.  The 
result  was  115  hits,  43  possible  hits  and  10  misses.  Or,  682  hits  and 
942  hits  or  possible  hits. 

C.  USING  SUBROUTINE  DET  FOR  CORRELATION 

The  same  set  of  test  signatures  used  in  VI-B  was  used  to  test  the 
programmed  correlation  and  selection.  The  input  and  results  are  listed 
in  Table  VI-3.  The  results  are  coded  as  follows: 
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Circuit 

Condition 

Run 

VV1 

1 

V1! 

Run  2 

VV1 

V1! 

R1 :  Short 

1 

1 

1 

1 

-90% 

1 

1 

1 

1 

-70% 

1 

1 

1 

1 

-50% 

1 

1 

1 

1 

-40% 

1 

1 

1 

1 

-20% 

1 

5 

3 

1 

+20% 

1 

3 

1 

3 

+30% 

1 

1 

1 

1 

+50% 

1 

2 

1 

1 

+60% 

1 

1 

1 

1 

+90% 

1 

1 

1 

1 

x3 

1 

1 

1 

1 

X75 

i 

1 

1 

1 

Open 

i 

1 

1  • 

* 

• 

LI :  Short 

i 

1 

1 

1 

-90% 

i 

1 

1 

1 

-70% 

l 

1 

2 

2 

-50% 

i 

1 

1 

1 

-40% 

i 

1 

1 

1 

-20% 

i 

1 

1 

1 

+20% 

l 

5 

1 

5 

+30% 

l 

3 

1 

1 

+50% 

l 

1 

1 

1 

+60% 

i 

1 

1 

1 

+90% 

l 

1 

1 

1 

X3 

i 

1 

1 

X75 

1 

1 

1 

Open 

l 

1 

1 

1 

R2:  Short 

3 

3 

3 

3 

-90% 

3 

5 

3 

2 

-70% 

3 

5 

3 

5 

-50% 

1 

1 

1 

1 

-40% 

1 

« 

1 

1 

1 

-20% 

1 

1 

1 

1 

+20% 

1 

n 

1 

1 

1 

+30% 

1 

1 

1 

1 

+50% 

1 

1 

1 

+60% 

1 

1 

1 

1 

+90% 

1 

1 

1 

1 

Open 

1 

1 

1 

1 

TABLE  VI-2. 

Results  Using  Sensitivity  Functions. 
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Cl :  Open  « 

-902  5 

-702  5 

-502  3 

-402  3 

-202  1 

+202  4 

+302  3 

+502  5 

+602  3 

+902  4 

X3  5 

X75  1 

Short  1 


i 

l 

l 

1 

1 

3 

5 

1 

1 

1 

1 

5 

1 

1 


«* 

i 

5 

5 

5 

3 

5 

5 

5 

5 

5 

3 

5 

1 

1 


* 

1 

1 

I 

1 

1 

3 
1 
1 
1 
1 
1 

4 
1 
1 


L2:  Short 
-902 
-702 
-502 
-402 
-202 
+202 
+302 
+502 
+602 
+902 
X3 
X75 
Open 


2 

1 

1 

1 

1 

2 

3 

5 

5 

5 

5 

5 

5 

1 


1 

5 

1 

4 
3 

5 
3 
3 

l 

i 

i 

i 

i 


2  1 

1  •  5 

1  5 

2  3 

1  3 

3  3 

4  3 

5  3 

3  1 

5  1 

5  1 

5  1 

5  1 

5  1 


C2:  Open  1 

-902  1 

-702  1 

-502  3 

-402  3 

-202  5 

+202  5 

+302  3 

+502  5 

+602  3 

+902  3 

X3  5 

X75  3 

Short  3 


1 

i 

• 

1 

1 

1 

3 

3 

1 

1 

1 

1 

1 

3 


1 

5 

5 

4 
3 

5 
5 
5 
5 
3 
5 
5 
3 
5 


1 

1 

1 

1 

3 

3 

3 

1 

1 

1 

» 

1 

3 

3 


TABLE  VI-2  (Continued) 
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1  -  correct  primary  selection 

2  -  one  of  two  primary  selections  count 

3  -  one  of  two  secondary  (no  primary  selections) 

4  -  incorrect  primary  selection. 

The  results  were  61  rated  1,  12  rated  2,  21  rated  3  and  2  rated  4. 
Or,  63%  completely  correct  and  98%  with  one  of  two  selections  correct. 
This  is  roughly  the  same  result  using  manual  correlation.  The  algorithm 
used  for  determining  primary  and  secondary  selection  is  outlined  in 
Section  VII-C. 


49 


Circuit  Condition 


VV1 


R1 

-90* 

-70% 

-60% 

+30% 

+60% 

+90% 

X3 

X75 


1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 


LI 

-90% 

-70% 

-60% 

+30% 

+60% 

+90% 

X3 

X75 


1 

1 

1 

2 

1 

1 

1 

1 


1 

1 

1 

2 

1 

1 

1 

1 


Cl 

-90%  3  3 

-70%  1  1 

-60%  3  3 

+30%  3  2 

+60%  3  2 

+90%  1  1 

X2  4  4 

X75  1  1 


R2 

-90%  1  1 

-70%  1  1 

-60%  1  1 

+30%  1  1 

+60%  1  1 

+90%  1  1 

X3  1  1 

X75  1  1 


TABLE  VI -3-  Computer  Correlation  Results. 


I  o 
k(. 

-90%  3  3 

- 70 %  3  1 

-60%  1  1 

+30%  2  2 

+60%  2  2 

+90%  3  3 

X3  3  3 

X75  3  3 


C2 

-90%  3  3 

-70%  2  2 

-60%  3  3 

+30%  1  1 

+60%  1  1 

+90%  1  1 

X3  3  3 

X75  2  2 


TABLE  VI -3  (continued) 


VII.  COMPUTER  PROGRAMS 


Several  computer  programs  ware  written  in  Fortran  IV  to  accomplish 
much  of  the  mathematical  work  associated  with  this  study.  Four  of  the 
programs  are  included  herein  with  detailed  descriptions.  Each  of  these 
programs  were  used  with  several  modifications.  However,  only  the  basic 
program  is  presented  in  each  instance.  All  listings  are  in  the  designated 
section. 

A.  PRIMARY  ANALYSIS  PROGRAM 

The  primary  analysis  program  contains  the  main  program  and  three 
subroutines.  A  fourth  subroutine,  PRQD,  is  used  for  root  finding  and 
is  not  included  in  the  listing. 

The  program  calculates  the  magnitude  of  both  network  functions  as  a 
function  of  the  complex  frequency,  s,  finds  the  roots  of  the  polynomials 
involved  and  selects  test  frequencies. 

Each  element  in  turn  is  then  set  to  each  out-of-tolerance  value  as 
selected  previously,  the  magnitude  of  each  function  at  each  test  fre¬ 
quency  is  calculated  and  all  responses  are  quantized  into  signatures. 

The  output  consorts  of  both  network  functions,  the  test  frequencies, 
the  nominal  response  at  each  test  frequency  and  the  list  of  library 
signatures. 

The  listing  is  the  program  which  calculates  the  test  frequencies  as 
explained  in  Section  IV-A.  To  modify  the  program  for  the  selection  from 
sensitivity  functions  the  portion  of  the  program  between  indicated  com¬ 
ment  cards  was  replaced  with  a  read  statement  and  the  test  frequencies 
were  used  as  the  input. 

Figure  2T  is  the  flow  graph  for  the  main  analysis  program. 
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Figure  21.  Flowgraph  of  Primary  Analysis  Program 
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B.  TWO  COMPONENT  VARIATION  PROGRAM 

The  main  analysis  program  was  modified  to  generate  signatures  for 
simultaneous  two-component  failure.  Two  different  modifications  were 
made,  one  generated  signatures  for  only  short  and  open  element  failures. 
The  other  allowed  for  short,  -50%,  +50%  and  open  failures. 

Since  the  programs  were  only  modifications  to  the  main  analysis 
program  flow  data  are  not  included  and  the  listings  are  omitted. 

C.  FAULTY  COMPONENT  SELECTION 

This  program  was  written  to  determine  possible  and  probable  component 
failures  from  the  two  input  signatures.  The  signature  from  V-j/I^  and 
Vg/V-j  are  both  processed.  A  component  is  flagged  as  possibly  failed  if 
one  t?  the  library  signatures  for  that  component  is  selected  as  a  match. 

If  two  or  mo^e  library  signatures  for  a  component  are  selected  as  a  match 
that  component  is  flagged  as  probably  failed. 

The  flow  graph  for  the  program  is  shown  in  Figures  22-24.  The  main 
program  allows  for  several  signatures  to  be  read  in  and  sent  to  subroutine 
DET,  one  pair  at  a  time. 

Subroutine  DET  first  checks  to  see  if  the  signature  indicates  the 
circuit  is  functioning  properly  (i.e.,  signature  is  all  5s).  If  so,  a 
message  is  printed  and  control  is  returned  to  the  main  program.  The  first 
entry  into  DET  with  a  non-nominal  signature  results  in  the  library  signa¬ 
tures  being  read  in  and  stored. 

When  it  has  been  determined  that  non-nominal  signatures  are  being 
processed,  both  input  signatures  are  compared  with  the  library  signatures 
to  determine  the  closest  match  distance.  The  signatures  are  compared  a 
second  time  and  those  library  signatures  that  are  the  minimum  distance 
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Figure  23.  Flowgraph  of  Subroutine  DET 
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front  the  input  are  used  to  set  flags  for  determining  possible  and 
probable  faulty  components  as  described  earlier. 

The  program  is  listed  showing  the  nearest  neighbor  matching 
procedure . 

D.  SUBROUTINE  LSQ 

Subroutine  LSQ  was  written  to  determine  which  library  signatures 
most  closely  matched  the  input  signature.  LSQ  is  similar  in  instruction 
flow  to  subroutine  DET  and  therefore  a  flow  chart  is  not  presented. 

In  order  to  determine  a  match,  LSQ  uses  a  method  similar  to  the 
nearest  neighbor  concept.  However,  if  one  considers  the  signature  to 
be  a  vector  with  each  test  frequency  a  dimension,  LSQ  does  not  weigh 
every  dimension  equally.  If  one  dimension  can  be  strongly  affected  by 
random  variations  in  component  values  that  dimension  is  not  considered 
equally  with  the  other  dimensions.  Thus,  the  distance  between  the  input 
vector  and  the  library  vectors  becomes: 

,!,  *i<*i  -  vs 

where 

t  =  b-jX-j  +  b2X2  +  b3X3  +  •  •  • 

=  input  vector 

Y  =  c-jYj  +  c2Y2  +  c3Y3  +  ••• 

=  library  vector 

a^  =  weighting  factor  for  the  i**1  dimension  (or 
frequency) 
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I.SQ  operates  on  only  one  signature  at  a  time  and  selects  the  proper 
library  signature  to  match  the  input.  The  output  is  a  list  of  the  pos¬ 
sible  signatures  and  the  component  variation  defined  by  each  one. 

One  may  choose  to  set  each  a^  =  1 .  This  would  be  the  case  if  the  ' 
quantisation  process  was  sufficient  to  take  care  of  the  random  variation 
prrbiem. 
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VIII.  SUGGESTED  FURTHER  STUDIES 

A.  ALTERNATE  APPROACH  TO  SIGNATURE  TABLE 

Table  VIII-1  lists  the  calculated  response  at  each  test  frequency 


for  the  primary  test  circuit. 

Table  YiII-2  lists  the  range  of  the 

response  at  each  frequency  and  the  maximum  deviation  from  nominal  response 

Zero  output  is  not  included. 

Test  Frequency  (Hz) 

Response 

VV1 

V1! 

.00048 

-2.?8  Db 

12.86  Ohms 

.027 

-6.77  Db 

3.22  Ohms 

.055 

-10.85  Db 

14.32  Ohms 

.299 

t56.8  Db 

3.01  Ohms 

2.99 

-120.3  Db 

4.34  Ohms 

49.97 

-215.5  Db 

53.46  Ohms 

TABLE  VIII-1 .  Nominal  Response  at  Test  Frequencies 


Test  Frequency  (Hz) 

Response  Range 

Maximum  Deviation 

v2/v1 (db) 

yx^o) 

V^V^db) 

yyn) 

.00048 

.0028 

*29.91 

3.0  66.77 

-27.63 

+53.91 

.027 

15.58 

-46.16 

.246  300.2 

-39.39 

+296.98 

.055 

-4.36 

-55.7 

3.0  302. 

-44.85 

+287,68 

.299 

-20.94 

-97.14 

,23  300. 

-40.34 

+296.99 

2.99 

-52.28 

•  -160.41 

3.  319.6 

+68.02 

+315.26 

49.97 

-98.48 

-255.5 

3.  5338. 

+117.02 

+5284.54 

TABLE  VIII-2.  Response  Range  and  Maximum  Deviation 

One  can  see  that  the  responses  do  not  cover  the  same  range  of  values  nor 
are  the  maximum  deviation  from  the  nominal  response  the  same.  From  the  nature 
of  the  network  functions  it  can  be  assumed  that  the  responses  at  different 
test  frequencies  are  not  independent  f  ;r  a  given  circuit  condition.  Due  to 

these  situations  it  would  seem  logical  to  weigh  equal  changes  at  different 
frequencies  differently. 
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The  setting  of  quantization  levels  and  subsequent  signature 
generation  is  one  way  to  accomplish  this  weighing.  The  setting  of 
quantization  levels  based  on  the  selected  situations  tends  to  weigh 
responses  at  some  test  frequencies  less  than  others.  However,  this 
Is  purely  an  "educated  guess"  method  and  has  no  analytical  foundation. 

An  alternate  scheme  was  studied  briefly  and  is  suggested  as  the 
subject  of  further  study. 

If  one  considers  each  test  frequency  as  one  of  the  dimensions  in 
an  n-dimensional  space  then  each  response  would  be  a  point  in  the  n- 
dimensicna!  space.  Instead  of  limiting  calculations  to  a  few  selected 
circuit  conditions  one  could  g'v.sr?te  many  points  representing  many 
circuit  conditions.  One  could  then  have  a  large  enough  sample  of  data 
to  dc  some  statistical  operations. 

In  order  to  remove  dependence  the  normal  procedure  is  to  multiply 
each  vector,  representing  a  point,  by  the  inverse  of  the  co-variance 
matrix.  Ihe  result  is  the  collection  of  points  in  a  "transformed  space" 
where  dependency  is  accounted  for. 

The  problem  is  the  construction  of  the  co-variance  matrix.  However, 
since  the  data  base  is  now  sufficiently  large, an  estimate  of  the  matrix 
can  be  calculated  from  the  data. 

With  a  circuit  under  test  the  response  would  constitute  a  vector 
which  would  be  multiplied  by  the  inverse  of  the  covariance  matrix.  At 
tMs  point  there  are  several  options  for  fault  identification.  One  can 
make  a  comparison  with  the  stored  data  and  make  the  fault  Identification 
using  the  nearest  neighbor  technique^  An  alternate  approach  would  be  to 
use  an  average  cf  the  nearest  m  neighbors  when  m  is  some  positive  integer. 
There  are  other  possible  procedures  for  matching. 
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If  one  is  only  interested  in  which  component  is  faulty  it  might 
also  be  of  interest  to  investigate  any  possible  grouping  of  points  which 
would  be  useful  as  in  pattern  recognition.  The  unknown  cotrtd  then  be. 
assigned  to  a  group  according  to  its  position.  Reference  5  details  the 
necessary  pattern  recognition  info  mat ion. 

Hie  above  described  procedure  would  require  a  computer  with  sufficient 
storage  to  process  and  store  the  information  pertaining  to  the  data  base 
and  covariatce  matrix.  However  such  a  computer  is  readily  available  for 
a  nominal  cost. 

B.  EXTENSION  OF  PROCEDURE 

It  soon  becomes  obvious  that  for  even  moderate  sized  circuits  the 
procedure  would  be  very  difficult  to  use.  For  example,  a  25-element 
circuit  with  6  possible  element  value  perturbations  would  yield  150 
signatures.  For  ease  in  testing  the  number  of  test  frequencies  should 
be  as  small  as  possible.  However,  with  150  signatures  considerable 
signature  duplication  is  probable  unless  the  number  of  test  frequencies 
Is  large. 

Reference  3  shows  how  a  ladder  network  may  be  divided  into  subnet¬ 
works  such  that  the  voltage  gain  of  a  given  subnetwork  is  independent  of 
the  voltage  gain  of  all  preceding  networks.  Therefore,  testing  by  add¬ 
ing  subnetworks  one  at  a  time  can  yield  the  fault  component. 

This  procedure  uses  only  the  voltage  gain  but  does  show  one  approach 
to  handirg  large  circuits.  In  any  case,  a  large  network  would  have  to 
be  divided  into  subnetworks  by  some  scheme  in  order  to  reduce  the  problem 
to  a  manageable  size. 
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IX.  CONCLUSIONS 


The  najority  of  the  conclusions  made  as  a  result  of  this  study  wer* 
based  on  a  comparison  of  the  results  is  section  VI  and  results  published 
by  other  investigators. 

Most  Investigators  reported  approximately  75X  isolation  using  one 
network  function  with  no  allowance  for  vari ati on- wi thi n-tol e ranee  of  the 
nonfaulty  elements.  One  investigator  [?]  reported  78X-98X  success.  How¬ 
ever,  the  results  were  not  sufficiently  documented  to  judge  the  merit  of 
the  claim. 

The  results  of  this  study  show  that  an  additional  error  rate  of 
172-42X  can  be  introduced  by  the  within  tolerance  variations  if  the 
conventional  procedures  are  followed.  However,  if  sensitivity  functions 
are  used  to  aid  in  selecting  test  frequencies  the  error  rate  due  to  other 
element  variations  drops  to  11X-31S  for  one  network  function.  Further¬ 
more,  using  two  network  functions  fault  isolation  is  successful  for  up 
to  S4X  of  random  single  fault  conditions. 

The  high  degree  of  successful  fault  isolation  clearly  points  out  the 
advisability  of  using  at  least  two  network  functions  for  testing.  The 
two  functions  used  for  this  study  were  Vg/V-j  and  V-j/Ij.  However,  the 
designed  purpose  of  the  circuit  may  dictate  a  different  choice  of 
functions  to  be  used. 

The  results  clearly  show  the  improved  accuracy  of  Isolation  using 
sensitivity  functions  as  a  guide  in  test  frequency  selection.  In  addition 
to  improved  isolation  this  method  also  eliminates  the  possibility  of 
unknowingly  picking  test  frequencies  which  yield  the  same  information. 
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Figure  26.  Signatures  using  Sensitivity  Procedure 
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Figure  27.  Signatures  for  Two  Component  Variations 
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COMPUTER  PROGRAM  LISTINGS 


REAL**  WMcf,A(320i 

REAL  KAG,N«AG,W«AG(  ?,320) 

niMcwcinM  riMcr.A  l 5 .8 » ,E( 1  O' ,D( « ) ,Qi 5) ,R ( 6  I .*( 5) ,UN( 51 . 
1YC  5  )  ,'F  I  ( *5 1  •■'»!  t  1  0  j  ♦  MAG (4,10;.  M-AG (6 , 1 0,  ?0)  ,PX (  6  )  , VX (  5 )  , 
22(5)  ,  X«AX(  ?0)»XMJM{20)  ,INMAf.f6,10,20).OMAG<2»105 
•  COMPCN  N^AG 

DATA  0MfGA/?E*0./, 0/^*0. /,R/5*0./, F 1/5*0. / ,F/5*0./ 

DATA  01/10*0./ 

READ  (5,110)  PAtO 
P{ 1 1=3. 

E( 2)=» l7 
E(3)=l  ,000001 
E(4>=10. 00001 
E{ 5 ) =9  . 

F( 6)=4. 

00  38  1=1,6 
00  38  K*l,10 

GO  TO  ( l, 2,?, 4, 5, 6, 7, 8, 9,10), K 

1  TEMP=E(  I) 

GO  TO  11 

2  F< I)=1.1*TPMP 
GO  TO  11 

3  E( I )=C.9*TEMP 
GO  TO  11 

4  F(I>=0. 

GO  TO  It 

5  E( 1 )=0.5*TFMp 
GO  TO  11 

6  E( I)=0.3*TEMP 

GO  TO  11 

7  E( I )  =  1  ,2*TP«P 
GO  TO  11 

9  E( 11=1  ,5*TFMP 
GO  T0  11 

9  E( I )=100.*TEMP 
GO  TO  11 

10  F( I )=TFVo 
cn  TO  39 

11  UN(1)=F(4) 

0(  1 )  =  F( 1 ) *F( 4 ) 

D(2)=F(2)*F(  5)+E(l  )*c(4)*(  P(3)  <E(6>) 

0(3)=F(5)*(T;(  1)*P(  3)*P(4)*P(6)  )+c(4)*E(2)*(F(3)^P(6H 
0(M=E<  5)*<E(l)*F(4)**rm*Ff6)+E(2l*F(3)) 

0( 5)=F( 2  )*E( 3)*F(4)*C(5)*E(6) 

Y(  1 ) =1  • 

Y(2)*P(4)*(P(3J*F(6)) 

Y<3)=E(3)*C(E) 

Y(4)=F(3)<E(4)*F(5)*C(6) 

IF  (K.GT. I.PR. T ,GT,1 )  GO  TO  36 
CC  WRITE  (8,100)  UN ( 1 ) 

CC  WRITS  (9,101) 

CC  WRITE  (8,102)  (0(N),N=1,5) 

CC  WRITE  (8,10?)  ( 0 ( N ) , N= 1,5) 

CC  WRITE  (P,10M 

CC  WRITE  (9,103)  (Y(N),M=1,4) 

C  — —  DELETE  STARTING  HFRE  FOR  SENSITIVITY  FUNCTION  USE 
DO  HU  LM  =  1 , 5 
DX  ( L  M | =0 ( L  M ) 

1111  YX(l«)=Y(LM) 

CALL  PROD  (0X,E,0,EI ,PCL,IR,IER)  • 

IF  (  1FR.C0,0)  GO  TO  12 
CC  WRITF  (8,105)  IEP 

CC  12  WRITE  (8,106)  (Q(N).EI(N) ,M=1,4) 

12  CONTINUE 

CALL  PROD  (YX,4,R,F,PnMtiptiER) 

IF  (  IER  ,c0  *0 )  GO  TO  13 
CC  WRITE  (9,105)  T E° 

CC  13  WRITE  (8,106)  ( P (N) ,F ( N) ,N=1 ,3 ) 

13  CONTI NUF 

DO  15  Msl,4 

? F  ( ABS(Pt(M)).LF. lc-10)  GO  TO  14 
OMEGA  <1,M)=ABS(EI(*) ) 
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GO  TO  15 

14  OMEGA  <  1  ,M)=A85I(Q(M)) 

»  »  »unr 

1  Z>  V^U^f  i  JM'fr 

00  l7  M=l,3 

IF  (  ABM  F(  m)  )  .LF.ie-10)  GO  TO  16 

OMEGA  (2,MJ=A8S(F<M)) 

GO  to  i? 

16  OMEGA  (2,M)-ABS(8(M|) 

17  CONTI MUE 

00  IB  v-  i » 2 
00  IB  M  =  1 » 4 
00  IB  NO  =  1 .4 

IF{OM=GA{MtN)8EO.OMEGA(MtNO)*ANO.N.NF.NO)  GO  TO  75 
GO  TO  13 

75  OMEGA  (m,N)=0. 

15  continue 

on  2i  j=t 1 4 

on  20  L  =  J « 4 

IF  fOMEGAfl*  JJ  .GT.OMFGAH  ,111  GO  TO  19 
GO  TO  20 

19  OMEGA  (1,5 )  =  O^EGA  (1,J) 

OMEGA  (1, JJ=OMcr,A  <1,1! 

OMEGA  (l,t)=OMFGA  (1*51 

20  CONT?  NtJF 

21  CONTINUE 

DO  22  1=1,4 

IF  ( OMEGA (  lfD.EO.0,1  GO  TO  22 

OMEGA  f  3,5>=<OM«-GA(l,lP)-OMEGAfl,t.n/2. 

OMEGA  (3,n  =  0MFGA  ( 1  ,L  )-0MFGA  ?3,5) 

OMEGA  (3,2)=QMFGA  ( t  ,i.  UOMcga  (3,5) 

IF  (0«EGA  (3,l).LT.O.)  nMFGA  ( 3 *  1 ) =OMEGA( 1  * U /2. 
LP2=L+2 

IF  (LP2,GT.4)  GO  TO  23 

ONEGA  (3t*l  =  (nMFGA(l,lP2>-OM=GA(l,LPn/2. 

OMEGA  (3,3)=OMEGA  (1 ♦ LP1 ♦OMEGA  (3,*) 

OMEGA  I 3,4)=0MCGA  ( 1 »Lp2 ) IOMEGA (3,5) 

IF  (LP2.EQ.4)  GO  TO  23 
LP3=t.  +  3 

OMEGA  (3,5)=(OMEGA  (1*IP3) -OMEGA  (l,LP2))/2. 

OMEGA  ( 3»4)*r'MFGA  ( 1  ,LP2 ) *0«EG4  (3,5) 

OMEGA  (3,5)=OMEGA(l,LP3)+OM6GA  (3,5) 

GO  TO  23 

22  CONTINUE 

23  00  26  J=1 ,3 
00  25  L=J»3 

IF  (OMEGA  (2,J).GT.0«EGA  (2,01  GO  Tp  24 
GO  TO  25 

24  OMEGA  (2,5)*0MEGA  (2,J) 

OMEGA  (2,J)=0MEGA  (2,U 

GMgGA  (2,LI=0MFGA  (2,5) 

25  CONTINUE 

26  C0NTINUC 
DO  2T  1*1,3 

IF  (OMEGA  (2,1) «E0.0. I  GO  TO  2? 

LP=L+1 

0**EGA  (4,5J  =  (0mEGA  (2,IP)-0«EGA  (2,U>/2. 

OMEGA  (4,l)=0MCGA  (2,l)-0MFGA  (4,5) 

.  OMEGA  (4,2 isOM^G A  (?,U+0«?0A  (4,5) 

OMEGA  (4,3)  =  OMFGA  (2,LP)+0MEGA  (4,*;) 

IF  IC«EGA(4,l?.LT.O, )  OMEGA  J 4 , l )  =  OMEG A ( 2 , L ) / 2 • 
LP2=L+2 

IF  UP2.GT.3)  GO  TP  23 

OMEGA  (4, 5 )= (0MFGA  ( 2,LP2)-0MEGA  (2*lP)|/2. 

OMFGA  ( 4, 3 )=OMEGA  t 2 ,l») +0MEGA  (4,5) 

OMEGA  (4,4)=0MtGA  ( 2 ,LP2 J+0«EGA  (4,5) 

GO  TP  28 

27  CONTINUE 

28  '2=0 
00  29  M=i,4 
00  29  Ml  =  3 ,4 
M2=M2*1 


NOT  REPRODUCIBLE 


*»e  miuii  ..nucrtiHi  .mi 

C  7  W  A  1  '<  I  *'  v*—»  '  A  *  I 

01 ( M2+ 1 ) =0MEGA( 3  *5 ) 

00  291  U  =  i,10 
00  291  IJ=IL,10 
IF  (OH  TJ).FQ.O.)  GO  TO  291 
ox*ni(!ii/oiu  J) 

IF  (OX.NF. 1.. AND. (OX. GT.. 9. AMD. OX. LT.l* 1)1 
291  CONTINUE 

00  32  M3=l,« 

00  31  V4=M3,a 

IF  (0i(M3).LF»01(N4) )  GO  Tn  31 
01 ( 10  *  =rl( M3 ) 

01(M3)=01( 

01  <M4j  ZOK  10) 

31  CONTINUE 

32  CONTI NUF 

00  34  M5=l ,9 

IF  !0l(M5).F0.0.l  GO  TO  34 
M6=l 0-M5 

IF  (M5.EQ.1)  GO  TO  35 
00  33  •*-»»!, 

33  01(M7)=01f M5+M7-1) 

GO  TO  35 

34  CONTINUE 

35  01(M6+U=0c 
M6=ha+i 

C  -  STOP  OFLFTION  HFPC 

36  CALL  MAGI  («6«ni ,UN,D,Y,M4G) 

IF  (K.NEd.OR.I.NE.l)  GO  TO  372 

00  37  LM=1,10 

MAG  (  3  ,LM)  =«AG(  l  ,L*) 

MAG  (4»LM)=MAG( 2«LM) 

37  WRITE  (4,115)  MS G { 3  » L M )  » MAG ( 4  ,  L  M ) 

CC  WRIT?  (8,10-')  (ni(MP),«3^1  ,mm 

WRITE  (4,115)  (01(M3),Mqil VM6) 

CC  WRITE  (S,10«) 

CC  WRIT?  ( P, 103 )  (MAG(t,N),MAG(2,N) ,N=1 ,M6) 

372  00  *00  !K«1,2 
DO  500  IL-1,«6 

IF(HAG(  IK.ID.LF.O.)  GO  TO  400 
DMAG(IK,IL)=20.*ALOG10(MAG(IK,ILI) 

GO  TO  500 
494  DMAG( IK,IL1=0. 

500  CONTINUE 

00  38  LM-l ,10 

IF  (MAG(3,LM).LF.O.)  GO  T0  371 

NMAG(I,K,IM)^»*AGU.IM)/MAG(3,IM) 

371  IF(MAG{4,LMHIF.Q.)  GO  Tn  *3 

NMAG( I ,K,LM+10)=MAG(2,LM) /MAG(4,LM> 

38  CONTINUE 

03  39  K=i,19 

39  CALL  NQM(6,2,3,K,XM4X(K) ,XMIN(K)) 

WRITE  (4,1151  (  X«AX(K7),XM!N(X.Z)  ,K?=1,19) 
00  40  K=l,9 

40  CALL  QUANTA  (K,XMIN(X)  ,XMAX(K),FACO) 

00  41  K=ii,19 

41  CALI  QUANTA  (K,XMIN(K)  ,XMAX(K),F4C0) 

00  411  1-1,6 

00  411  J=l«10 
00  411  K=i,?0 

411  INMAG(I,J,K)^NMAG(I,J,K) 

00  43  1=1,6 
00  4?  J=4*  9 

CC  42  WRITE  ( R,li2 )  ( INMAG{ I ,J,K) ,K=1,M6) 

42  CONTINUE 

CC  WRITE  (8,1131 

43  CONTINUc 

00  45  1=1,6 
DO  44  J=4» 9 
M6P=M6+10 

CC  44  WRITE  (8,112)  l!NMAG(I,J,Xi ,K=11 ,M6P) 

44  CONTINUE 


01 ( TL) =0* 
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;iR3r^ 


r  r 

uv 


uotte  Ifl.ini 
45  CONTI Ni'jF . 

WMTc  ( '*  ♦  1 1 2  )  I  (  I  INMAGI!  ,J,K1  ,K«ltM6)vJ:s4,9)»I«lv6) 
STOP  1 

100  F0RwAT  1 1 5X» 1PC1 3.5 ) 

101  FORMAT  <6X,*V2/V1  =  - 

l - * - M 

102  PORMAT  < 1 5X , 1PE 1 3 . 5 , 1 X ,*♦*, 1  PE  1 3 .5 , IX ,* $**1  +  *^13.5, 
HX,  •  c**?  ♦  • 

1,1PF.13.S1X,»S**3  ♦  »  ,li>E13.5,1  X,«S**4‘,/1 

103  format  t  I5x,tpri  3.5#*  +  •lDFi3„5,e  S**l  «■  ‘1PF13.5,*  S 
1**2  «•  *lPpT3.5*«  S**3 *  } 

104  F0PM.U  (4X,*  Il»l  - - 

1 - — — - — - *1 

10*  F0PMAT  U1GI 

106  FORMAT  (//// lOX.'PCiL  PART  I  »*AG  PART  REAL  PART  IMA*;  »A 
IRT  REAL  PART  I  MAG  PART  REAL  PART  t mag  PART* /OX.9F10.3J 

107  format  «////1x,'TK5  FOfCUFMClFS  USED  ARF:  * . /It>l0rl2.4) 

1  OR  *0RMAT  (//jx ,*V?/V!  MAGNITUDE  Zl»l  MAGNITUDE*  f\ 

109  «=0RM».T  (1PE1S.k,1PE21.5) 

110  FORMAT  <*5.0,151 

111  FORMAT  CtX,lPn20.5,lP2F20.51 

112  format  (1X,7T4! 

113  FORMAT  »1H  ) 

115  FORMAT  (2F15.7) 
lit,  FORMAT  (5E16.7) 

STOP 

6N0 


SUBROUTINE  MAGI  I N»0MEGA , XNUM* DEN,YNllMtMAG 1 
REAL  mag 

DIMENSION  OMEGA C 101 , XNUMI 5 1 »DEN(51 ,YNUM(5) ,MAGI4,101 
0?  1  1*1, N 
W=OMFGA< I ) 

YR=YNI«m»  \  J-YNUM(  3»  *W*W 
YI=YFilJM(2)*W-YNUMt41*W**3 
OR=DFN(  11-DENC  31  *M*W4-DEN( 5 1 *W**4 
D!=OENt  21*W-oc*.<4!*W**3 
YM*S0RT{ YP *YP+YI*YI 1 
OM=SORT(OR*nR+D!*DT ) 

NAG! 1,1 l=XNUMfl J/OM 
1  MAGI  2, 11=0*/ YM 
RETURN 
END 


SUBROUTINE  QUANTA  (K,OL,OH,FAC1 
C  K*TERM  IN  ARRAY 
C  OL* LOWER  LIMIT  FOR  TOl. 

C  OM=UPPER  C.IMIT  FOR  TOL. 

C  FACsFACTOR  FOR  QUANTIZATION 
REAL  NMAGI 6,10,20),ZI51 
COMMON  NM/-G 
DO  5  I  «l ,  6 
DO  5  J=ijlO 
QU^NMAGCt ,J,K) 

Z(li=PAC 

Z<2)=FAC*3. 

71 3)=FAC*6. 

ZI41=FAC*12. 

IF  I0U.LT.QL1  GO  TO  3 
Tc  IOU.GT.CH)  GO  TO  1 
NTAGI I • J»K)=0. 

GO  TO  5 

1  DO  2  L*l*4 

QFAC  =  ZIU*I0H-1.)i-Qm 
IF  tOU.GT.QFACl  GO  TO  2 
NMAG  I IfJ,K|=L 
GO  TO  5 

2  NMAG  I I,J,K)=5. 

GO  TO  5 
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•*  /A-  *  m  Vm  A  * »  *  ft*>-  Wxfl  »  " 


jjg"  J JJJJX I  uu.-^Ajg^y  u  I vuuuf 


f 


£ 

I 

k 


3 


4 

5 


00  4  M=1 «4 
0FAC=0L-Z<M)*< l. 

•  «*  «  At  •  t  »  Ar  % 

1  r  I  WUal  t  •  *‘Wf 


MU  A  n  f  f 


i/  \  .  UxA 


GO  TO  5  .  t/v 

NM&Gd ♦  J»K>«10. 

CONTINUE 

RETURN 

ENO 


-0t  ) 

r-rt  rn  &. 

*»w  *  ^  • 


SUBROUTINE  NO*  (  I  ♦  J1  •  J?  t  K *  X**AX  »  X*!  N) 
REAL  NMAG'^.10.20) 

COMMON  NMAG 
X“AX=1 . 

XM I N  =  1  . 

00  1  L=li! 

IF  (NKAG(LtMtK) .GT.XWJX) 

if  InmagilImNj.lt.xmin) 

CONTINUE 
R6TU°M 
END 


XM4X=NMAG(L t^t*) 
XM!N=NMAG(Lt^tKl 


1^  9  i^Kt.fH  •  A*.  •  Ak  •<•**«  «%  A  #  t  IM  •  A  t  ^  k 

tz  I  r'COiO  A  or\  JO{  J\;r.’7|0/  /  WO  1  O  l  O  I 

03  1  I«l/6 
09  1  <»l/9 

1  READ  <  5/ ICO )  (ISlGtI/K/L>#U*l/6)/USlG(I+6/K/U)/U»l#6) 
09  4  ! • 1/ 12 

09  4  J«  1  /  9 
09  3  <»  1  /  6 

3  JSI3(<)=ISI3( T/J/<) 

WRJTE(6#  1C1 )  !/J/(JSlG(Ki;/Kl«l/6) 

CAU,  USQ(  I#  J#  wS!3) 

4  CONTINUE 

100  F9fiVATC  1212) 

101  reRyAT(/#5Xi»THr  input  is:  ele«ent»/I3/ »  var.  nq»»/J2 
end 

SUBROUTINE  tS2< 11/ Jl/USIG) 

REAU  10  f  36) / I S 

DIMENSION  I  SI  3(36/6)/ JS1G(6)/KSI6(36/6)zFAC<6> 

IF  C 1 1-1)  1/6/1 

6  IF  (Jl.l)  1/9/1 
9  FAC(1)»1. 

FAC(2)«.3333 

FAC(2)»«3333<FAC(3)*»565jFAC(4)»*462jFAC!5)»*694 

FAC(6>*.745 

DO  c  J* 1/36 

2  READ {5, IOC)  (lSl3(J/I)/I*l/6)/(KS!3(J/I)/I»l#6) 

1  1F(IW)  3/7/7 

7  DP  5  U»l/36 
ID(J)»0 

09  4  1*1/6 

ID1*KSIG(J/I)-JSI3(;) 

4  ID<J)-ID<J)*FAC{I)*I01**2 

5  CONTINUE 
GO  TO  17 

3  D3  15  J*l/36 
IOC J) *0 

09  10  I»l/6 
I01«lSIG(J/I)*JSlG(n 
10  ID(J)*I0{J)*FAC<n»!Dl«*2 
15  CONTINUE 
17  IS«!D( 1 ) 

09  20  I «2/36 
1F(ID(J)-IS)  21/21/25 

21  lF(Il-7)  22/8/8 

8  WRITE(6/101)  U/ (<5IG( J/ I J)/ I J»l/6) 

GO  TO  25 

22  WMTE<6/1Q1>  J/(!SlG(U/lJ)/IJ?l/6) 

25  CONTINUE 

100  FORMAT (612/16/512) 

101  F8RMAT ( 5X/  *  A  P0SSI31.E  SIGNATURE  IS:  */ I3/5X/6I2) 

102  FORMATUH  ) 

WRITEI6/102) 

RETURN 

END 
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' '  kuruM,  >  i  ’»)W 


C  Ha  ini  PR  GGR  AM  ruR  USE  WITH  SUBROUTINE  DET  C  r  AC  NOT  INCLU 

DIMENSION  I$IG(12,<?,6)  ,  JS!6(  6)  ,LSIG(6) 

DO  l  1  =  1 » 6 
DO  1  K=1 ,0 

1  READ(2t  100)  <ISIGU,K,L),L=l,6),USIG(I<*6,K,L),L=l,6) 
DO  4  I=:,6 
DO  4  J=l,9 
DO  3  K=l,6 
JSIGIK)  =  !SIGU,J,K) 

3  LSIG<K)=IS!G{I+6,J,K) 

WRI  TE(3, 101 )  USIG(K)  ,K  =  1 , 6 )  *  {  LS IG(K )  ,K=1 ,6) 

4  CALL  DETU, J,JSIG,LSIG) 

STOP 

100  FORMAT C 1 2 1 2  I 

101  FORMAT ( /,5X- *  THE  INPUTS  ARE : • 14, 512 , • - *,612) 

END 


12 

3 


13 

23 

33 

1 

2 


16 

IT 

6 

? 


35 

36 

37 
10 


40 

45 

46 
15 

100 

102 

103 

101 


SUBROUTINE  CET(I,J,JSIG,LSIG) 

DIMENSION  ISIG(36,6),KSIG(36,6),JSIG(6),LSIG(6),ID(36) 
DO  3  K=i ,6 

IF(JSIG(K)-5)  13,12,13 
I FILSIGI K)-5)  13,3,13 
CONTINUE 
WRITEC  3 , 101 ) 

RETURN 

IF  C I  — 11  2,23,2 
IFU-l)  2,33,2 
DO  1  K=1 ,36 

RE  AD ( 2 , 1 00 )  (  ISIG(K,L),L=l,6i,(KSIGIK»L),L=l,6) 

DO  4  K=1 , 36 
ID(K)=C 
LD<K)=0 
DO  4  L-l*6 

ID1*= ISIGlK,L)-JSlGf  L) 

ID2=KSIG<K,L)-LSIG(L) 

I0(K)=ID(K)+101**2 
LD(K )=LD! K ) •MD2**2 
DO  5  L=1 ,6 
N I L ) =0 
IS=ID(1! 

ITsLOdl 
DO  7  11=2,36 
IFCLDC11I- 
IT=LD( III 
1FCIDU1I-IS) 

IS=ID(  ID 
CONTINUE 
DO  10  K=l«36 
IX*(K-l)/6+l 
IFUDCKI-IS)  35,35,36 
NCIX)=NI IXI+l 
IFILO(K)-IT)  37,37,10 
N(IX)=NUXD1 
CONTINUE 
DO  15  K=l»  6 
IF(NIK)-2)  40,45,45 
IFINIK)— 1)  15,46,15 
WRITEI3* 102 )  K,NtK) 

GO  TO  15 
WRITE! 3,103)  K 
CONTINUE 
RETURN 

F0RMATI6I2, 16,512) 


■IT)  16,17,17 


6,7,7 


FORM ATI* 
FORMAT  I  • 
FORM AT I « 
END 
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